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Abstract 

We discuss the kinetic theory of systems with long-range interactions. We contrast the mi- 
crocanonical description of an isolated Hamiltonian system described by the Liouville equation 
from the canonical description of a stochastically forced Brownian system described by the Fokker- 
Planck equation. We show that the mean-field approximation is exact in a proper thermodynamic 
limit. For N — * +oo, a Hamiltonian system is described by the Vlasov equation. In this collisionless 
regime, coherent structures can emerge from a process of violent relaxation. These metaequilibrium 
states, or quasi-stationary states (QSS), are stable stationary solutions of the Vlasov equation. To 
order 1/N, the collision term of a homogeneous system has the form of the Lenard-Balescu op- 
erator. It reduces to the Landau operator when collective effects are neglected. The statistical 
equilibrium state (Boltzmann) is obtained on a collisional timescale of order N or larger (when the 
Lenard-Balescu operator cancels out). We also consider the stochastic motion of a test particle in 
a bath of field particles and derive the general form of the Fokker-Planck equation describing the 
evolution of the velocity distribution of the test particle. The diffusion coefficient is anisotropic 
and depends on the velocity of the test particle. For Brownian systems, in the N — > +oo limit, 
the kinetic equation is a non-local Kramers equation. In the strong friction limit £ — > +oo, or for 
large times t S> £ , it reduces to a non-local Smoluchowski equation. We give explicit results for 
self-gravitating systems, two-dimensional vortices and for the HMF model. We also introduce a 
generalized class of stochastic processes and derive the corresponding generalized Fokker-Planck 
equations. We discuss how a notion of generalized thermodynamics can emerge in complex systems 
displaying anomalous diffusion. 
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I. INTRODUCTION 



The statistical mechanics of systems with long-range interactions is currently a topic of 
active research [1]. In a previous paper [2] (paper I), we have considered general models 
of Hamiltonian and Brownian systems with long-range interactions, for an arbitrary binary 
potential of interaction in D dimensions, and we have studied their equilibrium properties. 
In the present paper, we discuss the kinetic equations describing the out-of-equilibrium 
evolution of these systems. 

So far, most works have focused on the case of isolated Hamiltonian systems of particles 
in interaction such as the N-star problem in astrophysics [3, 4], the N- vortex problem in 2D 
hydrodynamics [4] and the HMF model [5, 6]. For these systems, the energy is conserved and 
the correct statistical ensemble is the microcanonical ensemble. In a proper thermodynamic 
limit iV — > +oo, the equilibrium one-body distribution function maximizes the Boltzmann 
entropy at fixed mass and energy. The "collisional" evolution of the distribution function 
is described by non-local kinetic equations such as the Landau-Poisson system (or the orbit 
averaged Fokker-Planck equation) in astrophysics [7] or the kinetic equations derived by 
Dubin & O'Neil [8] and Chavanis [9] for non-neutral plasmas and 2D point vortices. These 
equations increase the Boltzmann entropy Sb at fixed mass M (or circulation T) and energy 
E. For stellar systems, the relaxation time (Chandrasekhar time) scales as t re i ax ~ i^Tn^d, 
where is the dynamical time [7]. For the point vortex gas, the evolution is due to a 
condition of resonance which can be satisfied only if the profile of angular velocity is non- 
monotonic [8, 9]. When the profile of angular velocity is monotonic, the distribution is 
stationary on a timescale of at least Nt^ and it is not clear whether the system truly relaxes 
towards statistical equilibrium for longer times. The kinetic theory of the HMF model 
is also complicated [6, 10, 11] and seems to indicate a relaxation time scaling as t re i ax ~ 
iV L7 t£) [12]. In the "collisionless" regime, valid for sufficiently "short" times t <C t re i aX) 
the above-mentioned kinetic equations reduce to the Vlasov-Poisson and 2D Euler-Poisson 
systems. Since the relaxation time t re i ax increases rapidly (algebraically) with the number 
of particles, the Vlasov regime can be extremely long in practice (e.g., in astrophysics). The 
Vlasov-Poisson and 2D Euler-Poisson systems can undergo a phenomenon of collisionless 
relaxation (called violent relaxation in astrophysics) towards a metaequilibrium state on the 
coarse-grained scale [13, 14, 15, 16]. Since the relaxation time depends on N, the limits 
N — > +oo, t — > +oo (metaequilibrium) and t — > +oo, N — > +oo (statistical equilibrium) 
differ. The metaequilibrium state is described by non standard distribution functions which 
are nonlinearly dynamically stable stationary solutions of the Vlasov equation. In some 
cases, they maximize a H-function at fixed mass and energy [17, 18]. This maximization 
problem provides a refined criterion of nonlinear dynamical stability for the Vlasov-Poisson 
and 2D Euler-Poisson systems [6, 12, 19, 20, 21, 22, 23]. The statistical equilibrium state, 
reached for longer times, is described by the Boltzmann distribution. 

Since statistical ensembles are not equivalent for systems with long-range interactions, it 
is of interest, at a conceptual level, to introduce a canonical model of particles with long- 
range interactions. In that respect, we can consider a system of Brownian particles described 
by iV-coupled stochastic equations involving a friction and a random force in addition to 
long-range forces. These particles are in contact with a thermal bath that imposes the tem- 
perature T. For these systems, the correct statistical ensemble is the canonical ensemble. 
In a proper thermodynamic limit N — > +oo, the equilibrium one-body distribution function 
minimizes the Boltzmann free energy Fg = E — TSb at fixed mass M and temperature 
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T. The evolution of the distribution function of these Brownian particles is described by 
non-local Fokker-Planck equations which are the canonical counterpart of the kinetic equa- 
tions governing the evolution of Hamiltonian systems. For example, a gas of self-gravitating 
Brownian particles [24] is governed by the Kramers-Poisson and Smoluchowski-Poisson sys- 
tems. These equations decrease the Boltzmann free energy at fixed mass and temperature. 
Self-gravitating Brownian particles can experience an "isothermal collapse" [25], which is 
the canonical version of the "gravothermal catastrophe" [26] experienced by globular clus- 
ters. A Brownian model has also been introduced in the case of a cosinusoidal potential 
of interaction in d — 1 [6]. This is the canonical counterpart of the microcanonical HMF 
model. It could be called the BMF (Brownian Mean Field) model. 

In this paper, we compare these two descriptions (Hamiltonian and Brownian) and study 
their out-of-equilibrium properties. In Sec. II, we discuss the kinetic theory of Hamiltonian 
systems with long-range interactions by adapting the results of plasma physics to this more 
general context. In Sec. II A, starting from the Liouville equation, we consider a truncation 
of the BBGKY hierarchy in powers of the inverse particle number 1/N, which plays the same 
role as the plasma parameter in plasma physics. For iV — > +oo, the kinetic equation is the 
Vlasov equation. For long-range interactions, this equation is non-local due to mean-field 
effects and exhibits a "violent relaxation" towards a metaequilibrium state on a few dynam- 
ical times to (Sec. II B). In Sec. II C, we study the linear dynamical stability of a spatially 
homogeneous solution of the Vlasov equation and derive a criterion of stability generalizing 
the Jeans criterion in astrophysics. We also determine analytical expressions for the growth 
rate and damping rates of the perturbation. In Sec. II E, we study the nonlinear dynamical 
stability of a stationary solution of the Vlasov equation and compare with the Euler equation 
(Sec. II D). We also show that, for a spatially homogenous solution, the criterion of non- 
linear stability coincides with the criterion of linear stability. In the sequel of the paper, we 
consider a spatially homogeneous system which is stable with respect to the Vlasov equation 
and we study the time evolution of its velocity distribution function due to finite N effects 
(collisions). To order 1/N, the kinetic equation of a homogeneous system is the Landau 
equation (Sec. II F) when collective effects are ignored and the Lenard-Balescu equation 
(Sec. II G) when collective effects are properly accounted for. For 2D and 3D systems, these 
equations converge towards the statistical equilibrium state (Maxwellian distribution) on a 
relaxation time t re i ax ~ Ntu (for gravitational systems, the relaxation time is (N/ In N)tu 
due to logarithmic divergences). For one-dimensional systems, the Lenard-Balescu operator 
cancels out so that the relaxation, due to three-body (or higher) correlations, is longer than 
Nto (a similar result is obtained for the point vortex gas in two dimensions [8, 9]). In Sec. 
II H, we consider the relaxation of a test particle in a bath of field particles and derive the 
general form of the Fokker-Planck equation. The diffusion coefficient is anisotropic and de- 
pends on the velocity. This is responsible for anomalous diffusion and for a slow relaxation 
of the high velocity tail of the distribution [11, 27]. We provide various explicit expressions 
of the diffusion coefficient and friction force for a thermal bath with Maxwellian distribution 
function (subsections II HI and II H 2) and for one dimensional systems with an arbitrary 
distribution of the bath (subsection II H 3). In Sec. Ill, we study the temporal correlation 
function of the force and show that each mode decreases exponentially rapidly with a decay 
rate which coincides with the damping rate derived in the linear dynamical stability analysis 
of the Vlasov equation (Sec. II C). In Sec. II J, we consider the time evolution of the spatial 
correlation function of the particles in the linear regime and compare with the equilibrium 
results obtained in Paper I. In Sec. Ill, we develop the kinetic theory of Brownian systems 
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with long-range interactions. Starting from the N-body Fokker-Planck equation and using 
a mean-field approximation valid at the thermodynamic limit N — > +00, we derive a non- 
local Kramers equation. In the strong friction limit £ — > +00, or for large times t ^> £ _1 , 
it reduces to a non-local Smoluchowski equation. In Sec. Ill C, we study the evolution of 
the spatial correlation function for a Brownian system in the linear regime. Finally, in Sec. 
IV, we introduce a generalized class of stochastic processes and derive the corresponding 
generalized Fokker-Planck equations. We show that they display anomalous diffusion and 
that they are associated with a notion of generalized thermodynamics in /z-space. 

One interest of our general study is to present a unified description of systems with long- 
range interactions (Hamiltonian, Brownian, fluids,...) and to see how the results depend on 
the form of the potential of interaction and on the dimension of space d. Explicit results are 
given for gravitational systems, two-dimensional vortices and for the HMF model. Thus, 
our study shows the analogies and differences between these systems by placing them into 
a more general perspective. 



II. KINETIC THEORY OF HAMILTONIAN SYSTEMS 
A. The BBGKY hierarchy 

We wish to develop a kinetic theory of Hamiltonian systems with long-range interactions 
described by the iV-body equations (1-1) of paper I in order to obtain the evolution of the 
one-body distribution function f(r,v,t) = NmPi(r,v,t). We will see that many results of 
plasma physics developed for the Coulombian potential can be extended to a more general 
context. We shall discuss how these results depend on the dimension of space and on the 
form of the potential of interaction. We shall also discuss how the results are affected by 
the existence of a critical point in the case of attractive potentials (see Paper I). 

Starting from the Liouville equation (1-2), it is simple to construct the BBGKY hierarchy 
of equations for the reduced distribution functions, 
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f + E v-f + E E F(» ^HN-^J «ft* +I F + 1-0^ = 0. 

(1) 



i=l 1 i=l k=\,k^i i=l 



BP B f 

+ + (N - 1)— J F(2 - l)P 2 ( Xl , x 2 , t)d D * 2 = 0, (2) 



The first two equations of this hierarchy are 
BP X 

~~b7 

BP BP BP Bf 

~0f + vi ^ + F(2 - 1)^ + (N - 2)— J F(3 - 1)P3(X!, x 2 , x 3 , t)d D ^ + (1 - 2) = 0, 

(3) 

where x = (r, v). Introducing the decomposition (I-14)-(I-15), corresponding to the first 
terms of the Mayer expansion in plasma physics, they can be rewritten 

+iV— J F(2 ^ l)P^( Xl , x 2 , t)rf D x 2 = 0, (4) 
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+iV— y F(3 - l)P^(x 1 ,x 2 ,t)P 1 (x 3 ,t)rf D X3 
+ArA y F(3 - l)P^(x 2 ,X3,t)P 1 (x 1 ,t)rf D x 3 
- y F(3 -> l)P^(x!, x 2 , x 3 , t)rf D x 3 + (1 <-> 2) = 0, (5) 
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where we have taken N — 1 ~ N and N — 2 ~ iV for iV 1. We shall now consider the 
thermodynamic limit defined in Sec. II. B of paper I. For example, we can consider N — > +oo 
in such a way that the interaction potential (coupling constant) scales as ~ 1/iV", while 
f3 ~ 1, E /N ~ 1 and V ~ P 15 ~ 1. In that limit, the cumulant distribution functions 
P' scale as l/N j -\ We can therefore consider an expansion of the correlation functions in 
powers of the inverse particle number 1/N. This small parameter is the counterpart of the 
"plasma parameter" in plasma physics. Therefore, the methods of plasma physics can be 
applied in the present context with a different perspective. 



B. Vlasov equation and violent relaxation 

For N — > +oo, we get 

P 2 (x 1; x 2 ,t) = Pfa^Pifat) + 0(1/N) (6) 

so that the mean-field approximation is exact in a proper thermodynamic limit. In that 
case, the first equation of the BBGKY hierarchy reduces to 

df Of df 

W + V '^ +(F) -^ = ' (7) 
where we have introduced the mean-field force produced by the particles, 

(F)(r,0 = y F(l - 0)n( ri ,t)rf D ri . (8) 

The Vlasov equation (7) expresses the conservation of the distribution function in /i-space in 
the absence of "collisions" . In the case of long-range interactions, the Vlasov equation is cou- 
pled to Eq. (8). This coupling creates a complicated mixing process in phase space leading 
to a metaequilibrium state on a very short timescale, of the order of the dynamical time t D . 
This process is called violent relaxation in astrophysics [13] or chaotic mixing. It explains 
how a collisionless system can reach a quasi-equilibrium state (on a coarse-grained scale) 
as a result of phase mixing driven by long-range interactions. This metaequilibrium state, 
or quasi-stationary state (QSS), is a particular stationary solution of the Vlasov equation. 
Since it results from a turbulent mixing, it is particularly robust and has therefore nonlinear 
stability properties with respect to the collisionless dynamics. For a given initial condition, 
Lynden-Bell [13] has tried to predict the metaequilibrium state reached by the system by 
resorting to a new type of statistical mechanics taking into account the particularities of the 
collisionless dynamics (Casimirs) [28]. Unfortunately, his theory does not always give the 
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good result (i.e., the distribution that is actually reached by the system) because relaxation 
is incomplete in general so that the system can be trapped in a stationary solution of the 
Vlasov equation which is not the most mixed state. This concept of incomplete relaxation 
[13] explains why galaxies are more confined than predicted by Lynden-Bell's statistical me- 
chanics. In some cases, the metaequilibrium state turns out to maximize a H-function [17]: 

H[f] = - J C(f)d D rd D v, (9) 

where C is convex, at fixed mass M = J pd D r and energy E = \ J fv 2 d D rd D v + | J p$d D r. 
In the context of violent relaxation, Tsallis functional H q = —-^j J(f q — f)d D rd D v [29] 
is a particular H-function [18, 22, 23, 28], not an entropy, which occasionally, but not 
systematically, gives a good fit of the metaequilibrium state. Its maximization at fixed mass 
and energy leads to a particular class of nonlinearly dynamically stable stationary solutions 
of the Vlasov equation called stellar polytropes in astrophysics [23] . The same interpretation 
holds in 2D hydrodynamics [21, 30] and for the HMF model [6, 12]. We refer to [23, 28] for 
a more thorough discussion of these concepts. 



C. Linear dynamical stability 

We consider here the linear dynamical stability of a spatially homogeneous system with 
respect to the Vlasov equation (7)- (8). Following the standard derivation of the stability 
criterion and writing 5f ~ exp[i(k ■ r — cut)] with uj = u r + iX, the dispersion relation reads 
[6]: 

r k-^ 

e(k,w) = 1 - (2vr) D M(k) / ^d D v = (10) 

J k • v — ui 

where /(v) is the unperturbed distribution function and e(k, uS) is the dielectric function. 
The integration has to be performed by using the Landau contour [31]. Note first that in 
the case /(v) = p5(v), after an integration by parts, we obtain the explicit relation 

u 2 = (2n) D u(k)k 2 p. (11) 

On the other hand, for the Maxwellian distribution function with uniform density 

f(y)=(^Y /2 pe-^ i (12) 



2tt J 

the dielectric function can be expressed as 



e(k, u) = l + {2n) D u(k)(3mpW(J]M-) (13) 

k 



where 



2,7V 

is the W-function of plasma physics [31]. Explicitly, 



+ °° n, 2 



W(z) = -}= / e-^dx, (14) 



: y' .In 



W(z) = l — ze 2 / e^dy + iA —ze 2 . (15) 
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The dispersion relation can therefore be written 

1 + (2tt) d u(k)(3mpW(^/fM^) = 0. (16) 
We look for solutions of the form uj = iX where A is real. Then, we get 



where we have set r)(k) = — {2-n) D u(k)f3mp and defined the G-function 

G{X) = W^erfcM' (18) 

For x -> 0, G(x) = 1 + v/ttx + .... For x -> +oo, = 2x 2 (l + ^ + ...)■ For x -> -oo, 
~ — 27^ e-a:2 - We note that this function is always positive. Therefore, there exists 
solutions of the form uj = iX with A real only if u(k) < 0, i.e. for attracting potentials. 
For neutral plasmas with Maxwellian distribution, the pulsation u r of the perturbation can 
never vanish. 

The case of neutral stability uj = corresponds to r)(k) = 1. The case of instability 
(A > 0) corresponds to r)(k) > 1. In that case, the perturbation grows exponentially as 
5f ~ e xt . The criterion (17) can be written explicitly 



1 _ #) { 1 _^ erfc (^*)}.„, (19) 

where A is the growth rate. The case of stability (A < 0) corresponds to r](k) < 1. In that 
case, the perturbation decays exponentially as 5f ~ e~ 7 *. The damping rate 7 = —A is 
given by 

m-F(Ml) (20) 



2 k 



where we have defined the F-function 



F{X) l + v ^xe- 2 eT^ C (-x) , (21) 

such that F(x) = G(—x). For x — > 0, F(x) = 1 — y/nx + .... For x — > —00, F(x) = 
2x 2 (l + + •••)• F° r x ~^ +°°) F(x) ~ 2^. e ~ x ■ The criterion (20) can be written 
explicitly 



l-*){l + #^e rf o(-^2)}^0. (22) 

To summarize, the condition of linear dynamical stability is r)(k ) < 1 or, equivalently, 

1 + (2ii) D u(k)l3mp > 0. (23) 

It coincides with the criterion of thermodynamical stability (1-95) obtained in Paper I by 
considering the second order variations of entropy or the possible divergence of the spatial 
correlation function in Fourier space. The extension of these results to the case of polytropic 
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(Tsallis) distribution functions, instead of isothermal, is made in [6] for the HMF model and 
in [23] for self-gravitating systems. More generally, when / = f(v) depends only on the 
modulus of the velocity, we have df /dv = f'(v)v/v and the condition of linear dynamical 
stability is 

1 - {2n) D u{k) J ^d D v > 0. (24) 

The equality corresponds to the condition of marginal stability u = in Eq. (10). We refer 
to [6, 32] for a more thorough discussion of these results in the case of the HMF model. 



D. The Euler equation 

It is of interest to compare the stability of a system described by the Vlasov equation 
(7)- (8) to that of a system described by the Euler equations 

^ + V • (pu) = 0, (25) 



-^ + (u-V)u = — Vp-V$, (26) 
at p 

with a barotropic equation of state p = p(p) and $ = p * u. The energy functional of this 
barotropic gas is 

W[ P , u] = J p jT P -^dp'd D v + \ J p<S>d D r + f p^d D r. (27) 

This functional is conserved by the Euler equations. Considering the linear dynamical 
stability of a spatially homogeneous system with respect to the Euler equation, we find that 
the dispersion relation reads [6]: 

uj 2 = c 2 s k 2 + (2ix) D u(k)k 2 p, (28) 

where c 2 = dp/ dp is the velocity of sound. Note that for c s = 0, Eq. (28) coincides with Eq. 
(11). The condition of stability is 

c 2 + (2Ti) D u(k)p > 0. (29) 

In the unstable regime, u is imaginary (ou = iX with A > 0) and the perturbation grows 
exponentially as bp ~ e A *. In the stable regime, uj is real and the perturbation oscillates 
as 5p ~ e~ lwt . For an isothermal equation of state p = pT/m, the velocity of sound is 
c 2 = T/m and the instability criteria (23) and (29) coincide. In fact, this coincidence is 
general and goes beyond the isothermal distribution as shown in [6, 23] and in the next 
section. Indeed, for spatially homogeneous systems, the condition of dynamical stability 
(24) for a kinetic system described by the Vlasov equation with f — f(v) is the same as the 
condition of dynamical stability (29) for the corresponding barotropic gas (defined in Sec. 
II E) described by the Euler equation. Note, however, that the evolution of the perturbation 
is different in the two systems (kinetic and gaseous) [6]. 
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E. Nonlinear dynamical stability 



It can be shown that a distribution function /(r, v) which maximizes a H-function (9) at 
fixed mass and energy is a nonlinearly dynamically stable stationary solution of the Vlasov 
equation. This is because H, E and M are conserved by the Vlasov equation (see [23] for 
a more detailed discussion and references). Cancelling the first variations of H at fixed E 
and M, we get 

C'(f) = -(3e-a f = F((3e + a), (30) 

where F(x) = (C')^ 1 (—x). We note that / = /(e) depends only on the individual energy 
e = y + $(r) of the particles and is monotonically decreasing (/'( e ) < according to 
Eq. (30), assuming (5 > 0). This implies that the density p = f fd D v = p($) and the 
pressure p = J fv 2 d D \ = p(<E>) are functions of $. Eliminating $ between these two 
expressions, we find that the equation of state is barotropic in the sense that p = p(p). 
Therefore, to each kinetic system with distribution function / = /(e), we can associate 
a corresponding barotropic gas with the same equilibrium density distribution. Now, we 
have = f'(e)v 2 d D v = £ / §£ • vd D v = - / fd D v so that p'($) = -/>($) whi ch 

corresponds to the condition of hydrostatic equilibrium in a gas Vp = — pV$. Then, we 
get p'(p) = p'(<f>)/p'(<*>) = -p($ )/p'($ ). But, p'($) = / f'(e)d D v = j %/vd D v. Therefore, 
the velocity of sound c 2 s = p'(p) can be written c 2 s = —pj J ^/vd D v. This relation remains 
valid if the system is homogeneous [6] so that finally: 

c 2 = —^l . (31) 

J v 

Substituting this relation in Eq. (24) we see that the instability criteria (24) and (29) 
coincide as announced. 

We now turn to the nonlinear dynamical stability problem. It can be shown that a dis- 
tribution function /(r, v) which minimizes the Casimir-Energy functional F = E — TH 
(where T — 1//3 is a positive constant) at fixed mass is a nonlinearly dynamically stable 
stationary solution of the Vlasov equation. Again, this is because F and M are conserved 
by the Vlasov equation. If the distribution function minimizes F at fixed M (which is sim- 
ilar to a "canonical" stability criterion in thermodynamics) then it maximizes H at fixed E 
and M (which is similar to a "microcanonical" stability criterion in thermodynamics) [23]. 
However, the reciprocal is wrong in general so that the "canonical" criterion {min F \ M } 
is less refined than the "microcanonical" criterion {max H \ M,E }, and it just provides 
a sufficient condition of nonlinear dynamical stability. When the two criteria do not co- 
incide, this is similar to a situation of "ensemble inequivalence" in thermodynamics. Such 
"inequivalence" is observed for the nonlinear dynamical stability of self-gravitating systems 
such as stellar polytropes and is related to the Antonov first law [22, 23]. However, for spa- 
tially homogeneous systems, the two criteria are equivalent in general [6] and we shall use 
here the simpler "canonical" criterion. To minimize F[f] at fixed mass, we first minimize 
F[f] at fixed density profile p(r). This yields an optimal distribution /*(r, v) determined by 
C"(/*) = -P v i + A(r) where A(r) can be related to the density, using p = J f*d D v. Then, we 
minimize the functional F[p] = F[f*\ at fixed mass. It is shown in [22] that this functional 
can be written 

F[ P ] =Jpf P -jr d P' dD r + \J (32) 
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where p(p) is the equation of state of the corresponding barotropic gas defined above. As 
observed in [6, 22, 23], the functional (32) coincides with the energy functional (27) of a 
barotropic gas described by the Euler equations with u = 0. Now, a distribution p(r) 
which minimizes the energy functional (27) at fixed mass is a nonlinearly dynamically stable 
stationary solution of the Euler equations (25)-(26). This is because W[p, u] is conserved by 
the Euler equations. Since F[p] and W[p, 0] coincide, the condition of nonlinear dynamical 
stability for a homogeneous system described by the Vlasov equation is the same as the 
condition of nonlinear dynamical stability for the corresponding barotropic gas with respect 
to the Euler equation. This is the version of the nonlinear Antonov first law for homogeneous 
systems [6]. The second variations of Eq. (32) are 

6 2 p= f ^) ((5p) 2 rf D r + I [ 8p5M D^ (33) 

J 2p 2 J 

which must be positive for nonlinear dynamical stability. Now, repeating the same steps as 
in Sec. IV.D of Paper I by simply replacing T/m by c 2 s = p'(p), since p is constant, we find 
that the condition of nonlinear dynamical stability is 

c 2 s + (2n) D u(k)P > 0, (34) 
for all k. Using Eq. (31), this can also be expressed as 

1 - (2ir) D u(k) J f ~^-d D v > 0, (35) 

for all k. This generalizes the criterion of nonlinear dynamical stability given in [12] for 
the HMF model. Our approach also gives an interpretation of this criterion in terms of a 
condition on the velocity of sound in a barotropic gas with the same equation of state as 
the original kinetic system [6]. Indeed, we must have 

c 2 s > (c 2 s ) max = pv(k) max (36) 

for nonlinear dynamical stability. Finally, we note that the criteria (34) and (35) are the 
same as (29) and (24). This implies that, for homogeneous systems, the conditions of linear 
and nonlinear dynamical stability coincide. 



F. The Landau equation 

The Vlasov equation is very important for systems with long-range interactions because 
collisional effects become manifest on a timescale of order of Nt d or larger. For many realistic 
systems (e.g., galaxies in astrophysics) the number of particles is so large (N ~ 10 12 ) that 
only the collisionless regime matters for timescales of interest [7]. Therefore, we can consider 
the limit N — > +oo leading rigorously to the Vlasov equation. However, the limits N — > oo 
and t — > oo are not interchangeable. For sufficiently long times, collisional effects must be 
taken into account. This is the case for globular clusters in astrophysics [7] which form 
smaller groups of stars (N ~ 10 6 ) and whose age is of the order of the relaxation time. We 
would like now to take the effect of correlations between particles into account in order to 
describe the "collisional" relaxation due to finite N effects. In particular, we would like to 
obtain the form of the collision term to order 1/N. This is the first correction to the Vlasov 
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regime in the N^ 1 expansion of the correlation functions. It will describe the dynamics of 
the system on a timescale ~ Nto- 

There are different methods to obtain a kinetic equation for the distribution function 
f(r,v,t). One possibility is to close the BBGKY hierarchy by neglecting the cumulant of 
the three-body correlation function [31], which is of order P3 ~ N~ 2 in the present context. 
This corresponds to the Kirkwood approximation in plasma physics. Another possibility 
is to use the Klimontovich approach and develop a quasilinear theory (see, e.g., [3] and 
Appendix B). A third possibility is to use a projection operator formalism, e.g. [33]. An 
interest of this approach is that it takes into account non Markovian effects and spatial 
derealization. If we neglect collective effects, the projection operator formalism leads to a 
kinetic equation of the form 







x {^(1 - 0, t - t) A + ^(0 -+ 1, t - r) A J Z (ri , Vl , t - r)f(r, y,t-r). (37) 

Here, f(r,v,t) = NmPi(r,v,t) is the distribution function, (F)(r, t) is the (smooth) mean- 
field force and JF^(1 — > 0, t) = F^(l — > 0,t) — (F M )(r, t) is the fluctuating force created by 
particle 1 (with position and velocity r 1; Vi) on particle (with r, v) at time t. Between t and 
t — r, the particles are assumed to follow the trajectories determined by the slowly evolving 
mean-field (F)(r, t). Equation (37) is a non-Markovian integrodifferential equation. We 
insist on the fact that this equation is valid for an inhomogeneous system while the kinetic 
equations presented in the sequel will only apply to homogeneous systems. Unfortunately, 
Eq. (37) remains too complicated for practical purposes and we have to make simplifications. 
If we consider a spatially homogeneous system for which the distribution function / = /(v, t) 
depends only on the velocity [50], and if we implement a Markovian approximation, the 
foregoing equation reduces to 



df d 



dt dv» 



/+00 p 
dr J d D r 1 d D v 1 F^ l (1^0,t) 



x^d-o, t -,)(A_J!_)Z (Vl , i)/(v , t) , (38) 

with F M (1 — > 0,t) = F M (ri(t) — > r(t)). In our approximation, the particles follow linear 
trajectories with constant velocity since (F) = 0. Then, the collision term can be simplified 
(see Appendix A) and we obtain 

f = -( 2 -) D -7^ / d D v ld ^k^u(^ ■ (v - Vl) ] (4 - /||) . (39) 

where / = f(v,t) and f\ = f(vi,t). Introducing the relative velocity u = v — vi, this can 
be written more conveniently as 

d ! d w ( Jh^-M\ (40) 



dt dv» J v ; \ dv v dv\)' 

where 



K» v = 7r(2n) D m J d D kPA^(A;) 2 5(k • u). (41) 
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We note that the thermodynamic limit defined in Paper I amounts to considering that the 
coupling constant scales as 1/N, all other quantities being of order unity. Then, recalling 
that / = NmPi ~ N , we find that the collision operator in the Landau equation scales as 
1/N and represents therefore the first correction to the Vlasov equation. It also results from 
this observation that the Landau equation will describe the evolution of the system on a 
timescale of order Ntp. 

Equation (40) with Eq. (41) is the general form of the Landau equation. Landau derived 
it for an electronic plasma from the Boltzmann equation in a weak deflexion limit, using 
a linear trajectory approximation [34]. It can also be obtained from the Fokker-Planck 
equation by calculating the first and second moments of the velocity increments induced by 
a succession of two-body encounters [35] . The Landau equation can be further simplified by 
explicitly evaluating the tensor (41). In D = 3 and D = 2, we find that 

u^u v 



= K D ^ V -—j, (42) 
where 

i r+co i r+co 

K 3 = 87r 5 m- / k 3 u(k) 2 dk, K 2 = 87r 3 m- / k 2 u{k) 2 dk. (43) 

The Landau equation conserves the constants of motion of the Hamiltonian dynamics (mass, 
energy,...) and increases the Boltzmann entropy (H-theorem). For D = 2 and D = 3, the 
distribution function approaches the Maxwellian for t — > +oo. The relaxation time scales 
as t re i ax ~ Nto- Therefore, the kinetic theory justifies, at least at order 1/N and for ho- 
mogeneous systems, the Maxwell distribution predicted by the statistical theory. Since the 
microcanonical distribution is based on an assumption of ergodicity (see Paper I), it is not 
granted a priori that the system reaches statistical equilibrium. Thus, the development of a 
kinetic theory is necessary to validate (or not) the microcanonical distribution. For inhomo- 
geneous systems, it is not clear whether the kinetic equation (37) truly relaxes towards the 
mean-field Maxwell-Boltzmann distribution due to non-Markovian effects and spatial der- 
ealization. Therefore, the validity of the Boltzmann distribution depends on the dynamics 
and a kinetic theory is required to justify its relevance. 

The collisional evolution of stellar systems (that are inhomogeneous) is usually described 
by the Vlasov-Landau equation 

Of . .. Of ,™ Of f jD „ ml J, Of ,3h 



m ■ ■ it + (f) ■ ^ = j d \ ! ^- ! ^y m 

with now / = f(r,v,t) and f\ = f(r,Vi,t). This equation, obtained by combining Eqs. (7) 
and (40), assumes that the collisions can be treated as local (see [33] for a critical discussion 
of this approximation). It has to be coupled to the Poisson equation. For the gravitational 
potential in D = 3, 

K 3 = 2nmG 2 - / ^. (45) 
uj k 

This quantity exhibits a well-known logarithmic divergence at small and large scales. The 
divergence at small scales is common to both plasmas and gravitational systems. It is due 
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to the failure of the linear trajectory approximation. It is regularized by cutting the integral 
at the Landau length, corresponding to the impact parameter for large angle collisions. 
The divergence at large scales is specific to the gravitational case (in the plasma case, the 
integral must be cut-off at the Debye length as shown by the Lenard-Balescu treatment). 
It is due to the long-range nature of gravity and the absence of shielding. This problem 
has been the subject of several studies. In their stochastic analysis, Chandrasekhar & von 
Neumann [36] argue that the integral has to be cut-off at the interparticle distance since the 
distribution of the gravitational field (a Levy distribution) is dominated by the contribution 
of the nearest neighbor. Alternatively, most astrophysicists argue that the integral has to 
be cut-off at the size of the system, of the order of the Jeans length, which is presumably 
the gravitational analogue of the Debye length. In any case, because of the logarithmic 
divergence, the relaxation time scales as t re i ax ~ (N/ In N)t D instead of Nto- For the 
gravitational or Coulombian potential in D = 2, the integral behaves as J +o ° dk/k 2 and it 
diverges linearly for A = 2n/k — > +00. This suggests that K 2 is proportional to the large- 
scale cut-off (the Debye length in plasma physics). However, a more precise study based on 
the Lenard-Balescu equation taking into account collective effects is required to ascertain 
this result. 

In D = 1, the Landau equation becomes 



1 f d f 



J dvJ' 



where 



K = 4ir 2 m 



r+co 



ku(k) 2 dk. 



(46) 



(47) 



Due to the 5-function, the collision term vanishes identically. This implies that the distribu- 
tion function does not evolve on a timescale of order Nt D , i.e. df/dt = 0. Since the Landau 
equation is valid at order l/N, the relaxation towards statistical equilibrium will be due to 
non trivial correlations which appear at higher order in the large N expansion. This implies 
that, in ID, the relaxation time is longer than Ntu. For the HMF model, Yamaguchi et al. 
[12] report a non trivial scaling N 1 ' 7 . A similar cancellation of the collision term at order 
l/N occurs for 2D point vortices when the profile of angular velocity is monotonic. Indeed, 
neglecting collective effects, the kinetic equation describing the collisional evolution of an 
axisymmetric point vortex system is [9]: 



dP _N 1 2 d 
dt 4r dr 



r+00 

/ ndr^Q - fix) In 
Jo 



1-1^ 



r dr r\ dr\ 



(48) 



where P = P(r,t), P\ = P{r\,t) and r< (resp. r>) is the min (resp. max) of r and r\. 
The evolution is due to distant collisions. The collision operator is due to a resonance 
between vortices rotating with equal angular velocity Q(r) = il(ri) and vanishes identically 
when there is no resonance. Equation (48), derived in Chavanis [9] by using projection 
operator technics, is the vortex analogue of the Landau equation (40) for plasmas. A more 
general kinetic equation, taking into account collective effects, has been derived by Dubin & 
O'Neil [8] from the Klimontovich approach. It represents the analogue of the Lenard-Balescu 
equation (49) in plasma physics. As discussed above, the relaxation of 2D point vortices and 
of the HMF model toward statistical equilibrium is not clearly understood and demands to 
go to higher order in the l/N expansion. 
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G. The Lenard-Balescu equation 



The Lenard-Balescu equation can be obtained from the first two equations of the BBGKY 
hierarchy by neglecting non-trivial three-body correlations and assuming that the two- 
particle correlation function relaxes much faster than the one-particle distribution function 
(Bogoliubov's hypothesis) [31]. It can also be obtained from the Klimontovich equation by 
developing a quasilinear theory (see Appendix B). For a homogeneous system, the Lenard- 
Balescu equation can be written 

where 

r k- y~ 

e(k,u) = 1 + (27v) D u(k) / -^d D v, (50) 

j w-k-v 

is the dielectric function. The Lenard-Balescu equation can be seen as a generalization of 
the Landau equation taking into account collective effects. The classical Landau equation 
is recovered in the limit |e(k, k • v)| 2 ~ I. In D=l, the Lenard-Balescu equation becomes 



%-^h^&^-^-^) (51) 



where 



e(k,u) = 1 + 2nu(k) [ dv. (52) 

J uj/k — v 



uj/k — v 

We see that the collision term again vanishes identically at the order 1/N. 



H. Test particle in a thermal bath: the Fokker-Planck equation 

The Lenard-Balescu equation can also be used to describe the evolution of a test particle in 
a bath of field particles at equilibrium. In that case, we have to consider that the distribution 
fi of the bath is given, i.e. fi(v,t) = /o( v ) where /o( v ) is a stable stationary solution of 
the Vlasov equation (bath distribution). This procedure transforms the integro-differential 
equation (49) into a differential equation for the density probability P(v, t) of finding the 
test particle with velocity v at time t. It reads 

(53) 

where 

/k • ^ 
^d D v. (54) 
uj — k ■ v 
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Equation (53) can be written in the form of a Fokker-Planck equation 



dt dv^ \ dv u 

which involves a diffusion term 



" P " + (55, 



D"- = w(2n) D m J d D Vl d D WV ^ S{k ■ (v - vOl/ofa) (56) 



and a friction term 

\2 



rf = -n(2n) D m J .Pv/W f[k • (v - vi)]||(v 1 ). (57) 

The diffusion coefficient is due to the fluctuations of the force and is given by the Kubo 
formula 



r+oo 

D i»> = / (F^(0)F u (t))dt. (58) 
Jo 

The dynamical friction results from a polarization process and can be explicitly calculated 
by developing a linear response theory. In the present context, the coefficients of diffusion 
and friction depend on the velocity. Hence, it is more proper to write Eq. (55) in a form 
which is fully consistent with the general Fokker-Planck equation 

ep - 1 92 <p { ^)-±(rm. <») 



dt 2 dv*dv v V At J dv» \ At 
with 

At At dv v ' ' y J 

where the last equality is obtained from Eqs. (56) and (57) by using an integration by parts. 
We refer to Ichimaru [31] for a more comprehensive discussion of the test particle approach 
and for the connection between the Lenard-Balescu equation and the Fokker-Planck equation 
in the thermal bath approximation. Our aim, here, is to give explicit expressions of the 
Fokker-Planck equation and diffusion coefficient in particular cases. 



The isothermal case 



For the Boltzmann distribution 



MVi) = — 1 



D/2 2 

2.) 



(61) 



corresponding to statistical equilibrium (thermal bath), we easily find from Eqs. (56) and 
(57) that rf = (3mD^ u v u . Therefore, the friction coefficient is given by a generalized Einstein 
relation and the Fokker-Planck equation (55) takes the form 



dP 
~dt 



d 



(62) 
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This is similar to the Kramers equation in Brownian theory. However, in the present context, 
the diffusion coefficient D^ u is anisotropic and depends on the velocity v of the test particle. 
We note that for t — > +00, the velocity distribution of the test particle P(v,t) relaxes to 
the Maxwellian distribution (61) of the bath (thermalization). Writing 



5[k-(v- Vl )] = J e ik ' (v - Vl) ^, 



we can put the diffusion coefficient (56) in the form 



£>^(v) = (2tt) 



2D 



m 



+00 



dt I d D kk»k u 



u(kf 



_ik-vt 



/o(kt). 



(63) 



(64) 



J |e(k,k-v)| 2 

It will be shown in Sec. Ill that this equation can be interpreted as a Kubo formula (58) 
where t plays the role of time. For a Maxwellian distribution, the Fourier transform /o(kt) 
is a Gaussian. We can thus easily integrate on time t to obtain 



D^N) = 7r(27r) D mp( ^ 



1/2 



d D k 



kW m 2 o - P m 



k 



fk,k 



2k 2 



(65) 



On the other hand, for a Maxwellian distribution, the dielectric function is given by Eq. 
(13). Using (15), we can thus rewrite the Fokker-Planck equation in the form 



dP 
~dt 



1 d 

t R dx» 



G^(x 



f dP 
\dx v 



+ 2Px v 



where 



G^(x) = L D+1 



I dD] 



k^k^kf 



-(k-x) 2 



[l-rj(k)B(k 
and t R is a relaxation timescale defined by 



7T7/(A;) 2 (k-x) 2 e- 2 ( k - x ) 2 



t R l = 



(2n 



7T 

8D 



1/2 



m 



(Pm/2) 



1/2 



n L D+1 

v, k = k/k, B(x) 



(66) 



(67) 



(68) 



1 — 2xe x2 J* e y2 dy 



In the foregoing formulae, we have set x 
and r](k) = —(2ir) D u(k)Pmp. In addition, v m = (D/ ' (3m) l l 2 is the r.m.s. velocity and L is a 
lengthscale (size of the domain) introduced to make Eq. (67) dimensionless. The function 
B(x) can be written B(x) = 1 — 2xD(x) where D(x) = e~ x f*e y dy is Dawson's integral. 
We note the asymptotic behaviors B(x) = 1 — 2x 2 + ... for a: — > and B(x) 
x — > +00. Equations (66)-(68) provide the general form of the Fokker-Planck equation 
describing the evolution of a test particle in a thermal bath. They generalize the results 
obtained in [6, 10, 11] for the HMF model where the potential of interaction is truncated 
to one single mode and D — 1. We note that the relaxation time scales as t R ~ Nt D where 
we have introduced the dynamical time to ~ L/v m and the particle number N ~ nL D . As 
indicated previously, the Fokker-Planck collision term scales as 1/N in the thermodynamic 
limit defined in Paper I. 

We also note that in the case of the Coulombian interaction for which rj(k) = —k D /k 2 
where ko is the Debye wavenumber, the integral over k in Eq. (67) is of the form 



r+00 
Jo 



k D dk 



{B + k 2 /k 2 D ) 2 + C 2 



(69) 
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where B = £?(k-x) and C = y^lk-xje ( k ' x ) 2 are independent of k. The integral (69) can be 
rewritten I D = ^ +1 cVf D (5/C) with = t D dt/((x + t 2 ) 2 + 1). In the Landau 

approximation, we have instead J^ andau = / + °° k D ~ 4 dk which presents the divergences 
mentioned previously. In D = 3, the integral (69) is now convergent for k — > and, using the 
dominant approximation of Chandrasekhar [31], we have I3 ~ k%\n.{k m /kD) where k m is a 
small scale cut-off. Therefore, the Lenard-Balescu treatment justifies to cut the logarithmic 
divergence at large scales at the Debye length k^ 1 [31]. In D = 2, we note that the integral 
(69) is convergent for all values of k. The evaluation of the diffusion coefficient (67) in that 
case is left for a future study. Finally, the case D — 1 will be treated in Sec. II H 3. 



2. The Landau approximation 



If we make the Landau approximation |e(k, k • v)| 2 ~ 1 which amounts to neglecting 
collective effects, we can rewrite the Fokker-Planck equation in the form 



dP 1 d 



where now 



dt t R dx» 

G""(x) = / rf D k^^ e - (k - x)2 , 



° V m JO 



(70) 



(71) 



(72) 



Equation (71) also represents the asymptotic behavior of the diffusion coefficient (67) for 
|v| — > +00 and for T — > +00 since in these limits e ~ 1. In that case, the tensor G^ v 
given by Eq. (71) can be calculated explicitly by introducing spherical or polar systems of 
coordinates. In D = 3, 



x^x 



CT = (G^--G ± )- S - + -G ± 6^, 



Xr 



with 



2vr 3 / 2 
G11 = G(x), 



x 



2tt 3 / 2 

G\ = [etf(x) - G(x)\, 



x 



where 



G(x) = 



2 1 



1 



ert(x) -=e 

v 71 " 



(73) 



(74) 



(75) 



2 f x _ 2 
erf(x) = —= / e * dt. 

V 71 " Jo 



(76) 



In D = 2, 



CP'" = (6*|| -G ± )^ + G ± 5^, 



(77) 
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with 



ire 



G 



ire 



(78) 



In the above expressions G\\ and Gj_ are the diffusion coefficients in the directions parallel 
and perpendicular to the velocity of the test particle. Finally, in D = 1, 



G(x) = 2e 



(79) 



We can also develop a test particle approach in the case of 2D point vortices to describes 
the stochastic evolution of a test vortex in a bath of field vortices at statistical equilibrium 
[4, 9, 37]. Making a thermal bath approximation, we can transform the kinetic equation 
(48) into a Fokker-Planck equation of the form 



dP 
~dt 



ld_ 

r dr 



dP 

rD(r)(— + (3 1 P ^ 



dlpeq 



with a diffusion coefficient 



8 Sr 



\nN(u) eq (r), 



(80) 



(81) 



rQ! (r) is the local shear created by the statistical distribution of field vortices 



where E(r, 

(the angular velocity fi eg (r) is related to the vorticity by uj eq {r) = (l/r)(r 2 Q eq )'). In Eq. 
(80), P(r,t) denotes the density probability of finding the test vortex in r at time t. In 
the more general case where the field vortices have a distribution wo(r) which is a stable 
stationary solution of the 2D Euler equation (not necessarily the equilibrium Boltzmann 
distribution), the Fokker-Planck equation (80) is replaced by 



dP 
~dt 



ld_ 

r dr 



n( \f dP v dlnuj o 



(82) 



and the diffusion coefficient is still given by Eq. (81) with u Q in place of (cu) eq . In particular, 
for a vorticity profile u>o(r) = Ae~ Xr of the field vortices, it is easy to see that the diffusion 
coefficient of the test vortex decreases like D(r) ~ r 2 g -Ar 2 £ Qr r _^ _|_ 00 _ yy e emphasize 
the analogies with the Fokker-Planck equations derived previously in the case of material 
particles. Note that for point vortices, the dynamical friction is replaced by a systematic 
drift [37] along the vorticity gradient. 



3. The case D = 1 



In D = 1, the Fokker-Planck equation (53) can be written 



OP _ d 
dt dv 



dP rfln/o 
D(v)[ — P— — 

ov dv 



where D(v) is given by 



D(v) = An 2 mU(v) / dk 



ku{k) 2 
\e(k, kv)\' 



(83) 



(84) 



18 



Equation (83) is similar to a Fokker-Planck equation describing the motion of a Brownian 
particle in a potential U(v) = — ln/o(v) created by the field particles. In D — 1, the 
distribution function of the test particle P(v , t) always relaxes toward the distribution of 
the bath fo(v) for t — > +00 (the typical timescale governing the approach of the test particle 
to the bath distribution is ~ Ntc). In other dimensions, we see from Eq. (53) that this 
is true only when /o( v ) is the Maxwellian distribution (statistical equilibrium state). The 
particularity of the dimension D — 1 is that the Lenard-Balescu collision operator (49) 
cancels out for any distribution function f(v) while in D = 2 and D = 3 the cancellation 
of the collision operator occurs only when the distribution is Maxwellian. In D = 2 and 
D = 3, an arbitrary distribution of the field particles /o(v, £) changes on a timescale of 
order Ntr> as it relaxes to the Maxwellian due to the development of correlations (finite N 
effects). Therefore, the test particle approach is valid only for t <C Ntc, in an interval of 
times where we can consider that / (v,t) is approximately stationary. But in that case, the 
distribution of the test particle P(v,t) does not relax toward /o(v). It is only when /o( v ) 
is the Maxwellian that we can consider that the distribution of the field particles is frozen 
(for any time). In that case, P(v,t) relaxes toward /o(v) i n a time ~ Nt^. The situation is 
different in D — 1. In D — 1, a stable stationary solution of the Vlasov equation (any) does 
not change on a timescale of order Nto- Therefore, we can consider that the distribution of 
the field particles fo{v) is frozen on a time ~ Nto (or larger) which is precisely the time it 
takes to a test particle to relax toward fo{v). The fact that the distribution P(v, t) of a test 
particle relaxes to any distribution fo(v) of the field particles explains why the distribution 
of the field particles does not change on a time ~ Nt D . 
For the Maxwellian distribution in D — 1, one has 



dP _ 1 d 

dt t R dx 



fdP 
G W (- + 2Px 



(85) 



with 

G(x) =2L 2 [ °° dk- — kT]( ; k / e X — — - T , (86) 

These results extend those obtained in [6, 10, 11] for the HMF model for which 

2L 2 e~ x2 

G{X) = (T/T c -B(x)) 2 + nx 2 e- 2 * 2 - (88) 

Note that for periodic potentials, the integral over k must be replaced by a discrete 
summation over the different modes. On the other hand, for a Coulombian potential 
r/(k) = —k%/k 2 , the integration on k can be made explicitly, using &i(x) = n/A — tan _1 (a;)/2 
and we obtain 



2 x 



1 tan 



71 \^/7^\x\e 



(89) 



with asymptotic behaviors G(x) ~ ^/TiL 2 k 2 D /\x\ for x — > ±00 and G(0) = L 2 A;|). 
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I. Temporal correlation functions 



We now wish to determine the temporal correlation function of the force experienced by 
the test particle. If we ignore collective effects, the temporal correlation function is given by 

(F"(0)F"(t)) = N J F»(l -> 0, 0)F"(1 -> 0, t)P 1 (r 1 , Vi)d D r lC 2 D Vi. (90) 

Using Fourier transforms and making a linear trajectory approximation as in Appendix A, 
we obtain 

(F"(0)F"(i)) = m(2vr) D J Wk v u<& 2 e^™f Q (\ 1 )d D \ 1 d D \i. (91) 
We can also write Eq. (91) in the form 

(F"(0)F^)) = m(27r) 2D J PA^kfe"^ f (kt)d D k. (92) 
For a Maxwellian distribution, we get 

/2 2 
A^A^k^e-^VWcPk. (93) 

For the gravitational interaction, we can easily perform the integrations by introducing a 
spherical system of coordinates. The correlation function can finally be written 

C^ = (C ll -±C ± )^f+ 1 -C ± 5^, (94) 

with 

AnpmG 2 

11 = — Vt — ^' ^ ' 

C± = 47rp ™° 2 [erf ( x ) - G(x)}. (96) 
where x = (/3m/2) 1 / 2 v and the function G(x) is defined by Eq. (75). We note the result 

(F(0)-F(t)) = ^^erf(x), (97) 

which shows that the correlation function of the gravitational force decreases as i -1 . This 
asymptotic result was first noted by Chandrasekhar [38] using a different approach. The 
diffusion coefficient is then obtained from the Kubo formula (58) by integrating over time 
t. This returns Eqs. (73)-(75) except that the logarithmic divergence J dk/k in Eq. (72) 
now appears on the time integration f dt/t. The divergence at large scales is connected to 
the very slow (algebraic) decay of the temporal correlation function. This algebraic decay 
is strikingly in contrast with usual Markov processes in which the correlations decrease 
exponentially rapidly with time. Coming back to Eq. (93), we note that each mode has a 
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gaussian decay ~ exp[— k 2 t 2 /2f3m] but the integration over all modes leads to an algebraic 
behavior ~ £ _1 . 

If we come back to the general problem and take into account collective effects, it is 
shown in Appendix B that 



(F"(0)F"(*)) = m(2n) D J 



u(kf 



ikvd 



f (^)d D Vl d D k. 



e(k, k • vi) 



(98) 



If we integrate over time t and use the Kubo formula (58), we recover the expression (56) 
of the diffusion coefficient. Introducing 



Q(k,t) 



e^ , ,/° (vi) ,J \ 



|6(k,k. Vl )| 2 

we can write the correlation function more compactly as 



(F"(0)F"(i)) = m(2n) D I k"k v u(k) 2 e- ik - vt Q(k, t)d u k. 



Vf, /1,\2„— ik-vt / 



(99) 



(100) 



For a Maxwellian distribution (61), 



'Bm\ l,2 \ r +(X1 p- pm %? 



(101) 



where we have chosen k in the direction of x and set uj = kv x . To determine the large time 
asymptotics of Eq. (101), we need to investigate the poles of the function 



e 2fe 2 



\e\ 2 (k } u) 

Using the notations introduced in Eq. (67), we have 

1 2 



(102) 



\e\ 2 (k,iu) 



l-rj(k)B 



(3m uj 



IT 



UJ 



\k 



(3muj 



(103) 



Setting uj = iX where A is real, we find after straightforward algebra that the foregoing 
expression can be rewritten 



\e\ 2 (k,iX) = e(k,iX)e(k, —iX), 



where we recall that 



e(k,iX) = 1 - r)(k)/G 



(3m X 
~Yk 



(104) 



(105) 



Equation (102) can thus be rewritten 



e(k, iX)e(k, —iX) 



(106) 



21 



Clearly, this is an even function of A. We need to determine the values of A for which the 
denominator vanishes. Since the distribution of the bath is stable by hypothesis (T > T c ), 
we have r](k) < 1 for all k. Therefore, the denominator cancels out for A = ±7 where 
7 > is determined by e(k, — 27) = 0. Thus 7 coincides with the damping rate of the stable 
perturbed solutions of the Vlasov equation (see Sec. II C). It is given as a function of k by 
Eq. (20). Next, we consider A = ±7 + A where A C 1. Expanding Eq. (106) for small 
values of A, we find after elementary calculations that for u — > ±27 



where the constant is given by 



2 h 2 1 

K(l) = ~ (108) 

Since f(u) behaves like a Lorentzian for uj — > ±ry, this implies that, for t — > +00, the 
different modes in the correlation function decrease exponentially rapidly as 

Q(k, t) ~ e"^, (109) 

V P"^7fc |_p/( J§ELJk^ 

with a rate 

7* = (HO) 

depending on the wave vector A; (see Eq. (20)). These results generalize those obtained in 
[6, 10] for the HMF model. In this case, the correlation function decreases like 

(^(0)^(0) - "8^ Vln^)l C ° s( ^ )e ' (1U) 

where the decay rate 

7 = (2/f3) 1 / 2 F'\ V ) (112) 

only depends on the temperature (recall that 77 = /3kM/4ir). For T — > T c + , 7 ~ 
(8/A;M) 1 / 2 (T - T c ) and for T -> +00, 7 ~ V2TlnT (to leading order). Close to the critical 
temperature, the correlation function decreases very slowly. iVote that this slow decay may 
invalidate the Markovian approximation close to the critical point and lead to dynamical 
anomalies. On the other hand, at high temperatures, the decay is very fast. In fact, had we 
ignored collective effects and used Eq. (93), we would have obtained a gaussian decay 

(F(0)F(t)) = ^-co8(vt)e~^, (113) 
An 

instead of an exponential decay with large 7-rate. This shows that collective effects are 
important even far from the critical point since they modify the large time behavior of the 
temporal correlation function. For the gravitational potential, Eq. (110) becomes 



Wi*Wf )■ (H4) 



(3m \ k 2 
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Using F l (x) = ^=(1 — x) + ... for x — > 1 and F 1 (x) ~ V - In a; for x — > 0, we find that 



2 - / *3 



t^V^H 1 -^' (115) 

7fc ~__ln(- , (A;^+oo). (116) 



\[Wn \kj 

J. Evolution of the spatial correlation function in the linear regime 

To conclude this section on the kinetic theory of Hamiltonian systems with long-range 
interactions, we would like to discuss the time evolution of the spatial correlation function. 
Returning to the second equation of the BBGKY hierarchy (5) for the two-body correlation 
function and considering the case of a homogeneous medium, we obtain 

BPL BP' .dP 1 , , 

_^ +Vl _^ +F(2 ^ 1) P l(V2) _i (Vl ) 

+N ^I F ^^ 1 'l P ^^t)Pi^i)d D ^ + (l"2)^0, (117) 

while the first equation (4) gives BP\/Bt = for — > +oo. We assume that, initially, no 
correlation is present among the particles. Then, for sufficiently small times (linear regime), 
the correlations will be small and we can neglect the integrals in Eq. (117). This yields 

a pi a pi f)p 

-gf + vx^ + F(2 - 1)^^)^^) + (1 <- 2) = 0. (118) 

We now assume that, initially, Pi ~ exp[— j3mv 2 /2] is the Maxwellian distribution. As we 
have seen previously, this distribution is conserved to leading order at later times. Intro- 
ducing the correlation function h though the defining relation 

P 2 (ri,vi,r 2 ,v 2 ,f) =Pi(vi)Pi(v 2 )[l + /i(r 1 -r 2 ,Vi -v 2 ,t)], (119) 

we find that it satisfies the differential equation 

^ + u-^ = /5mF(2^1)-u, (120) 

where x = ri — r 2 and u = vi — v 2 . For a stationary solution, we recover Eq. (1-53) of 
Paper I. Note also that Eq. (1-51) can be obtained from the static (B/Bt = 0) expression of 
Eq. (117) by keeping the integrals and assuming that Pi is Maxwellian. 

Returning to Eq. (120), taking its Fourier transform and solving the resulting first order 
differential equation, we get 

h(k, u, t) = -I3m 2 u{k){l - e- lk - ut ), (121) 

where we have assumed that, initially, the system is uncorrelated. Taking the inverse Fourier 
transform of Eq. (121), we finally obtain 

/i(x,u,f) = /3m 2 [w(x- uf) -u(x)]. (122) 
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It may be useful to discuss the HMF model explicitly. In that case, Eq. (120) becomes 



dh dh 

m +u d4> 



k/3 

u sin i 

2tt 



The stationary solution of Eq. (123) is 



K4>) 



kp 
2k 



COS0. 



(123) 



(124) 



This is similar to Eq. (1-66) of Paper I, but it does not display the expected divergence as 
we approach the critical point T c . This is due to the neglect of the integrals in Eq. (117). 
Furthermore, it is noteworthy that the solution of Eq. (123) does not converge towards the 
stationary distribution (124) for t — > +oo. Indeed, the explicit time-dependent solution of 
Eq. (123) is 



, k(3 
h((p,u,t) = — 



cos 4> — cos(0 — ut) 



(125) 



which shows an oscillatory behavior. This is not really surprising since Eq. (125) is only 
valid for short times. In this linear regime 



kd 

h((f>, u, t) = — — sin(0) ut (t — > 0). 



(126) 



III. KINETIC THEORY OF BROWNIAN SYSTEMS 
A. The non-local Kramers equation 

We shall now derive the kinetic equations describing a system of Brownian particles with 
long-range interactions defined by the stochastic equations (1-39)- (1-40) of Paper I. We shall 
use and generalize the method of Martzel & Aslangul [39, 40]. We start from the general 
Markov process 

P N (r 1 ,v 1 ,...,r N ,v N ,t + At) = J d D (Ar 1 )d D (Av 1 )...d D (Ar N )d D (Av N ) 

xP iV (r 1 - Ari, vi - Avi, ...,r N - Ar N ,v N - Av N ,t) 
xw(t! - Ari, vi - Avi, - Ar^, - Avjv|Ari, Av x , Ar^, Av N ). (127) 

where w denotes the transition probability from one state to the other specified by the term 
in parenthesis. Using Eq. (1-39), it can be rewritten 

w(r 1 , vi, ...,r N , v N \ Ari, Avi, ...,Ar N , Av N ) = 
5(An - viAi)...5(Arjv - v JV A*)^(ri, vi, r N , v N \ Avi, A Vjv ). (128) 

Then, the integration over Ari...ArAr is straightforward and yields 

P N (T 1 ,v 1 ,...,r N ,v N ,t + At) = J rf D (A Vl )...rf D (A Vjv ) 

xP N (r l - viAt, vi - Avi, ...,rjv - v N At, v N - Av N ,t) 
x^(ri - viAi, vi - Avi, r N - v N At, v N - Av N \ Avi, Avtv), (129) 
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or equivalently 

P JV (r 1 + ViA*,Vi,...,r JV + v JV A*,v JV ,i + A<) = J d D {Aw 1 )...d D {Aw N ) 

xPjv(ri, vi - Avi, ...,rjv, vjv - Av N ,t) 
x^(ri, vi - Avi, r N , v w - Avjv| Avi, Avjy). (130) 

Expanding the right hand side in Taylor series and introducing the Kramers-Moyal moments 

1 



M, 



ni...njv 



>i, vi, r N , vat) = lim Ai -o- 



- j d D (Av 1 )...d D (Av N ) 



At ni\...riN' 

x (-A Vl ) ni ...(-A ViV ) nAr V(ri, vi, r N , v N \ A Vl , A Vjv ), 



we get 



N 



dP N ^ <9P, 



A' 



E 



9" 



Qn N 



n N 
N 



{M ni _„ nN P N }, 



(131) 



(132) 



dt ^ dri ^ dv^'"dv' 

i=l ni...njv 

where the sum runs over all indices such that X^n^ — 1- For the stochastic process (1-39)- 
(1-40), only a few moments do not vanish, namely 



Af ...ni=i...o = + mV i C/(r 1 , r^), 



(133) 



Mr 



0...n»=2...0 



P. 



(134) 



Substituting these results in Eq. (132), we obtain the iV-body Fokker-Planck equation 

N , OT , or-, \ JV 



dP N dP N dP N 

z=i ^ 



9* 



<9ri 



f9v, 



e- 



i=l 



(135) 



where Fj = — mVjf/(r 1 , r^) is the force by unit of mass acting on the i-ih particle. The 
TV-body Fokker-Planck equation decreases the free energy 



F[P N ] = (E)[P N ]-TS[P N ] 



(136) 



constructed with the average energy (1-8) and the entropy (1-26) defined in Paper I. Indeed, 
an explicit calculation yields 

F = -j^j ^^D d ^ + iP N ^\ D v 1 d D ^ l ..A D v N d D ^ N < 0. (137) 

Since F = at equilibrium, the term in bracket in Eq. (135) vanishes by virtue of Eq. 
(137). Since d/dt = 0, the advective term must also vanish, independently. From these two 
requirements, we find that the stationary solution of the iV-body Fokker-Planck equation 
corresponds to the canonical distribution (1-42) which minimizes the free energy. 

We can now obtain the equivalent of the BBGKY hierarchy for the reduced distribution 
functions. It reads 



dt 



t^ + t t nk^^ + (N-j)±[d^ +1 F( 3 + l 

i=l 1 i=l k=l,k^i 1 i=l J 







dPi 



dvi 



E- 



dP 

^ + 6^(138) 
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In particular, the first equation of the hierarchy is 



dt 



+ Vl |i + (J v_ 1) / W ^ 1) ^ = ^ 



avi 



(139) 



We can similarly obtain an equation for the two-body correlation function. Implementing the 
decomposition (I-14)-(I-15) of Paper I, neglecting the cumulant of the three-body correlation 
function, and taking the thermodynamical limit defined in Paper I, we find that P' 2 ~ 1/N. 
Therefore, in the limit N — > +00, we can make the mean-field approximation 



P 2 (ri, vi,r 2 , v 2 ,t) = Pi(ri, Vi,*)Pi(r 2 , v 2 ,t). 



We thus obtain 



where 



dPi dP x dP x 



d 



dt 



dri 



dvi 



(F) 



-Nm 



J d D r 2 d D w 2 



dvi 

du 



dvi 

ri - r 2 )Pi(r 2 , v 2 ,f) 



Introducing the distribution function / = NmPi, this can be rewritten 

d_ 
dv 



df df 9/ 
— + v • — + (F) • — 



where 



(F) = -V$ = - J d D r> P(r',t)^(r-r') 



(140) 
(141) 
(142) 

(143) 
(144) 



is the mean-field force acting on a particle and we have set £ = Df3m. Equation (143) is a 
non-local Kramers equation. It decreases the free energy 



F[f]=E[f]-TS[f] 



fv 2 d D rd D v + ^- I p^d u r + T I f\nfd n rd u v, 



(145) 



which plays the role of a Lyapunov functional [18]. The free energy (145) is the result- 
ing expression of Eq. (136) in the mean-field approximation. We note that the non-local 
Kramers equation (143) is obtained at the same level of approximation (i.e. for iV — > +00) 
as the Vlasov equation (7) for Hamiltonian systems. The introduction of a friction and a 
random force in the equations of motion (1-39)- (1-40) yields a "collision term" of the Fokker- 
Planck form in the right hand side of Eq. (143). This "collision" term selects the mean- field 
Maxwell-Boltzmann equilibrium distribution (1-24) among the infinite class of stationary 
solutions of the Vlasov equation (left hand side). This mean- field Maxwell-Boltzmann dis- 
tribution extremizes the free energy (145) at fixed mass. Furthermore, only minima of free 
energy are linearly dynamically stable via Eq. (143) [18]. We should also contrast the non- 
local Fokker-Planck equation (143) to the local Fokker-Planck Eq. (62). Although they look 
similar, their physical content is quite different. Indeed, Eq. (55) describes the motion of 
a single test particle in a thermal bath at statistical equilibrium while Eq. (143) describes 
the evolution of the whole system of N Brownian particles in interaction out of equilibrium. 
Therefore, Eq. (55) is a local differential equation while Eq. (143) is non-local due to the 
mean-field force produced by the distribution function f(r,v,t) evolving in time. In addi- 
tion, in the present case, the diffusion coefficient D is given in Eq. (1-40) while in the test 
particle approach it is derived from the Hamiltonian dynamics at the order 1/N for TV » 1. 
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B. The non-local Smoluchowski equation 

In the strong friction limit £ — > +00, or equivalently for large times t ^> it is possible 
to neglect the inertial term in Eq. (1-40). In that case, we are led to consider a system of N 
Brownian particles in interaction described by the coupled stochastic equations in physical 
space 

^1 = - f xm 2 V i U(r 1 ,...,r N ) + y/2D*Ri(t), (146) 
at 

where \x = l/m£ is the mobility and D* = D/£ 2 = T/m^ is the diffusion coefficient in 
physical space. The Einstein relation reads \x = (3D*. We can now repeat the above 
procedure to obtain the kinetic equations governing the evolution of the Brownian system 
in the overdamped regime. For a general Markovian process in physical space, we have 

P N ( ri , ...,r N ,t + At) = J d D (Ar 1 )...d D (Ar N ) P N (n - An, r N - Ar N , t) 

xw(ri - An, r N - Arjv| Ari, Ar N ) (147) 

where w denotes the transition probability. Expanding the right hand side in Taylor series 
and introducing the Kramers-Moyal moments 



-^ni...njv 



we get 



(r 1 ,...,r^) = lim A ^ -— / d D (Ar l )...d D (Ar N ) (-A^.^-Ar^ 

At ni\...n N \ J 

xw(r h rjv| Ari, ArM}148) 



dP N d ni d nN 



E ^■■■»3»{Mn 1 ...n N PN} (149) 



dt ^ drl 1 Or 



■ N 

«i...njv 



where the sum runs over all indices such that J2i n i ^ 1- For the stochastic process (146), 
only a few moments do not vanish, namely 

M ... ni= i...o = fim 2 V l U{r 1 ,...,r N ), (150) 



■\ A,..,, •>..,, = A, (151) 
Substituting these results in Eq. (149), we obtain the iV-body Fokker-Planck equation 



dR 



N 



N 



dt 



i=i 



_d_ 



8Pn ,n 9 TT , . 



(152) 



The stationary solutions of this equation correspond to the configurational part of the canon- 
ical distribution (1-44). Again, it is possible to derive the equivalent of the BBGKY hierarchy 
for the reduced distribution functions. It reads 



dPj _ A d 



i=i 



dr. 



+ \xm 



E Pj^ + >un%X-j) 



k=l,k^=i 



dri 



P, 



du 



i+i" 



'^ +1 d D r 



j'+i 



. (153) 
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Introducing the decomposition (I-14)-(I-15), neglecting the three-body correlation function 
and considering only terms of order 1/N or larger in the thermodynamic limit N — > +00, 
we obtain 



dP l 
~dt 



dt 



_d_ 
dri 

_d_ 



D^ + vLNrrtP^t) J P 1 (r 2 ,t)^d D r 2 + pNm 2 J P^n,r 2 ,t)^d D r 2 

d P' du r 

^*^r + ^ m2p i( r i^) p i( r 2^) -^r + Nfim 2 P^n,r 2 ,t) I Pi(r 3 , 



+Npm 2 Pi(r 1 ,t) J P'(r 2 ,r 3 ,t)^d D r 3 



(154) 

t) d ~^d°r 3 
dri 

+ (1 <-> 2)55) 



The stationary solutions of these equations coincide with the equations of the equilibrium 
BBGKY-like hierarchy (see Paper I) as expected. In the limit N — > 00, we can make the 
mean-field approximation 

P 2 (n, r 2 , t) = Pi(n, t)Pi(r 2 , t). (156) 
The equation for the density then reduces to 

dPi _d_ 
dri 



which can be written 



dt 

dp 
dt 



n dPi n <9$ 

D* — h pmP 1 — 

dri dri 



(157) 



where $(r, t) is related to p(r,t) as in Eq. (144) 
equation. This equation decreases the free energy 

F[P\ = \ 



V- [D*(Vp + pmpV <!>)}, (158) 
This is the non-local Smoluchowski 



p$d D r + T I p\npd u r 



(159) 



which plays the role of a Lyapunov functional. The Smoluchowski equation can also be 
derived directly from the Kramers equation in the strong friction limit £ — > +00 [41]. This 
can be done by working out the moments equations of Eq. (143) as in [18] or by using a 
Chapman- Enskog expansion as in [42]. To leading order in l/£, the velocity distribution is 
Maxwellian 

'Bm\ D/2 

- p(r,i)e-^ +0(r i), (160) 



and the evolution of p{r,t) is given by Eq. (158). The free energy (159) can be obtained 
from Eq. (145) by using the approximate expression (160) of the distribution function to 
express the free energy as a functional of p. 

Considering the linear dynamical stability of a spatially homogeneous distribution of 
particles with respect to the non-local Smoluchowski equation (158) with $ = p * u, we 
immediately get the dispersion relation 

i£u = —k 2 + (2ii) D u(k)k 2 p. (161) 
m 

We note that the condition of instability (corresponding to ioo < 0) is the same as for an 
isothermal distribution described by the Euler or the Vlasov equation (see Sees. II C and 
II D). However, the evolution of the perturbation is different. In the unstable regime uj = iX 
with A > 0, the perturbation grows exponentially rapidly as 5p ~ e xt . In the stable regime 



CO 



—27 with 7 > 0, the perturbation decreases exponentially rapidly as 8p 



~ e 



-it 
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C. Evolution of the spatial correlations 



In this section, we study the development of the spatial correlations (at the order 1/N) 
for a Brownian system in the overdamped regime. For a homogeneous system, Eq. (155) 
for the two-body correlation function reduces to 



dh _d_ 
dt dr 1 

or, equivalently, 

dh 
dt 



dh 2 dui2 

D *~x h V m ~a — 

or i or i 



+ pmp I h(r 2 -r 3 ,t)-^-d D r 3 



■I 



+ (l<-2), 



2£LA 



h(x,t) + f3m 2 u(x) + (3mp J h(y, t)u(x - y)d D y 



(162) 



(163) 



(164) 



(165) 



where h eq {k) is the equilibrium value of the correlation function (in Fourier space) given by 
Eq. (1-54) of Paper I. For the BMF model (one mode), we find that 



where x = r x — r 2 . Taking its Fourier transform, we get 
dh 



+ 2D^[l + {2ti) d (3mpu{k)]h = -2D*(3m 2 k 2 u(k) 



dt 



This equation is easily integrated in time to yield 

h(k,t) = h eq (k)\l - e -2^ 2 [l+(2^/W«( fc )]t 



h(<f),t) = 



(3k/2n 



1 _ e -2D*(l-/3//3 c )t 



COS0. 



(166) 



For t — > +oo, the correlation function relaxes towards its equilibrium form (1-66). Note 
however that the relaxation time diverges as we approach the critical point since t re i ax ~ 

(1-/V&)- 1 . 



IV. GENERALIZED KINETIC EQUATIONS AND EFFECTIVE THERMODY- 
NAMICS 

A. Generalized Kramers and Smoluchowski equations 

We shall introduce a class of stochastic processes leading to generalized Fokker-Planck 
equations. These equations are associated with an effective thermodynamical formalism in 
/x-space. Let us first consider the case of non-interacting Langevin particles described by 
the stochastic process 

dr 

— = -AiV$e*t(r) + y/K(T,t)R(t), (167) 

where R(t) is a white noise and Q ext (r) an external potential. Since the function in front of 
R(t) depends on the position, the last term in Eq. (167) can be interpreted as a multiplicative 
noise. The corresponding Fokker-Planck equation is 

^ = ±A(K(r,t)p) +/iV(pV$ e ^). (168) 
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When K is constant, one recovers the usual Fokker-Planck equation with a diffusion coeffi- 
cient D = K/2. Recently, there was some interest for the nonlinear Fokker-Planck equation 

^ = A(Dpi) + pV(pV$ ext ), (f 69) 

that arises in connection with Tsallis generalized thermodynamics, see e.g. [43]. Defining 
a generalized temperature T = 1/(3 through the Einstein-like relation T = D/p, we can 
rewrite the foregoing equation in the form 

^ = V[D(Vp ? + PM]. (170) 

It is easy to verify that the nonlinear Fokker-Planck equation decreases the Lyapunov func- 
tional 

F[p\ = J p$ ext d D v + ^J (jfl - p)d D r, (171) 

which can be interpreted as a free energy F = E — TS associated with the Tsallis entropy 
S q = — ^-j- / (p q — p)d D r. Furthermore, the stationary solutions of this equation are given 
by a g-distribution 



P 



(172) 



Q 

which minimizes the Tsallis free energy at fixed mass. In an attempt to justify the nonlinear 
Fokker-Planck equation from a microscopic model, Borland [44] proposed to consider the 
generalized stochastic process 

- = -//V<t> e:ct (r) + V2Dp{v, t)<*-WBJ(t). (173) 

The last term is a multiplicative noise which depends on r and t through the density p(r, t). 
Therefore, there is a feedback from the macroscopic dynamics. For this stochastic process, 
K(r,t) = 2Dp q ~ 1 and the Fokker-Planck equation (168) takes the form (169) where the 
diffusion coefficient depends on the density. 

In a previous work [18], we remarked that the nonlinear Fokker-Planck equation (170) 
is a particular case of a larger class of generalized Fokker-Planck equations associated with 
generalized entropy functionals encompassing Tsallis entropy. Similar observations have 
been made independently by Frank [45] and Kaniadakis [46]. These equations can be written 

^ = V{D[pC"(p)V P + P P V^ ext }}, (174) 
where C is a convex function, and they monotonically decrease the functional 

F = J p<$> ext d D v + T J C(p)d D r, (175) 

which can be interpreted as a generalized free energy. When C(p) = plnp, we recover 
the usual Fokker-Planck equation associated with the Boltzmann free energy and when 
C(p) = (p q — p)/(q — 1) we recover the nonlinear Fokker-Planck equation associated with 
the Tsallis free energy. Comparing Eq. (174) with Eq. (168) we find that 

p = DP, \v{Kp) = DpC"(p)Vp. (176) 
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The second equation is equivalent to 



Integrating by parts, we get 



l -Kp = J P DpC"(p)dp. 



(177) 



-Kp = D[pC'(p)-C(p)], 



(178) 



so that, finally, 



K(r,t) = 2Dp 



C{p) 



(179) 



Therefore, a stochastic process leading to the generalized Fokker-Planck equation (174) is 
given by 



dr 

— = _^v$«rt(r) + \\2Dp 



C{p) 



R(t), 



(180) 



where R(£) is a white noise. We note in particular that the strength of the stochastic force 
depends on the local density of particles. When C(p) = plnp, we recover the usual Langevin 
equation and when C(p) = (p q — p)/(q — 1) we recover the stochastic process considered by 
Borland [44]. The same arguments can be developed in velocity space (taking into account 
the inertia of the particles), see [47], leading to the generalized Kramers equation 



dt dr ext <9v dv 



D[fC"(f)^ + Pfv 



181) 



B. Generalized non-local Kramers and Smoluchowski equations 



We shall now consider the generalization of the previous approach to the case of particles 
in interaction. We thus consider the N coupled stochastic equations 



dTj 
dt 



Vi, 



(182) 



dt 



-£vj - VjC/(ri, r N ) + J 2Dfi 



C(fi 
fi 



Ri(t), 



(183) 



where fi = /(rj, Vj, t). We shall say that these equations describe a gas of Langevin particles 
in interaction [48]. We reserve the term of Brownian particles for usual diffusion when 
C(f) = / In / is the Boltzmann function (see Sec. Ill A). Using exactly the same steps as 
in Sec. Ill A, we can derive the generalized iV-body Fokker-Planck equation 



N 



dP N ^ dP N 

1=1 v 



dt 



dR 



N 



dri 



N a 

Y — 



C(fi 
fi 



Pn + ZP N v t 



(184) 
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Implementing a mean-field approximation, we arrive at the generalized non-local Kramers 
equation 



df df ,^,df d 
— + v— + (F)— = — 
dt dr dv dv 



D[fC"(f)^ + (3fv 



(185) 



where (F) is the mean-field force (144) and we have set £ = Df3. This equation decreases 
the Lyapunov functional 



F[f\ = \ I fv 2 d D vd D v+ l - I P *d D r 1 I C(f)d D rd L \. 



(186) 



which can be interpreted as a generalized free energy. 

Considering the strong friction limit £ — > +oo, or the limit of large times t ^> we 
can derive a generalized Smoluchowski equation. This is done in [18] from the moment 
equations of Eq. (185) and in [42] by using a Chapman-Enskog expansion. The generalized 
Smoluchowski equation can be written 



dt 



-(Vp + pV$) 



(187) 



The fluid is barotropic in the sense that p = p(p) where the equation of state is entirely speci- 
fied by the function C(f). For the Boltzmann entropy, we recover the ordinary Smoluchowski 
equation (158) with p = p-^. The generalized Smoluchowski equation (187) decreases the 
Lyapunov functional 



F[ P \ 



(188) 



which can be interpreted as a generalized free energy. It can be derived from the generalized 
free energy (186) by using the leading order expression of the velocity distribution C'(f) = 
— /3[^-+A(r, i)]+0(l/£) to express F[f] as a functional of p (see [42] for details). Considering 
the dynamical stability of a spatially homogeneous solution of the generalized Smoluchowski 
equation (187) with $ = p * u, we get the dispersion relation 



iuj£ = elk 2 + {2n) D u{k)k 2 p, 



(189) 



with c 2 s = p'(p) which generalizes (161). For the damped Euler equations [18], the term iuj^ 
is replaced by uo{uo + it;). For £ = we recover (28) and for £ — > +oo we recover (189). 

We can also obtain a form of generalized Smoluchowski equation by starting directly from 
the coupled stochastic equations in physical space 



drj 
dt 



-pV i U(r 1 , ...,r N ) + \ 2D*pi 



C(A) 



Ri(t), 



(190) 



where pi = p{Yi,t) and R»(t) is a white noise acting independently on each particle. We 
note that these equations cannot be obtained from Eqs. (182)-(183) by simply neglecting 
the inertial term, as is done in the standard Brownian case. Applying steps similar to those 
of Sec. Ill B, we obtain the generalized iV-body Fokker-Planck equation 



dPj 



N 



dt 



N 



V 



d 



dri 



d_ 
dri 



D*pi 



C( Pi 



N 



d 

+ pP N ^U(r 1 , ...,r N ) 



(191) 
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Implementing a mean-field approximation, we finally get the generalized non-local Smolu- 
chowski equation 

^ = V[D.(pC"(p)Vp + PpV*)], (192) 

where <£> is related to p as in Eq. (144) and we have set p = D*f3. This equation decreases 
the Lyapunov functional 

F[P] = \f P®d D r + T J C(p)d D r, (193) 

which can be interpreted as an effective free energy. 

We note that the stationary solution of the generalized Fokker-Planck equation (191) is 
determined by the integro-differential equation 

A d 

P N \ + pPn—U^, r N ) = 0, (194) 

where p{ji) = f P N d D r 1 ...d D r i _ 1 d D r i+ i...d D r N . Contrary to the case of Brownian particles 
studied in Sec. Ill B, the equilibrium iV-body distribution does not seem to have a simple 
form. In particular, it does not seem to be possible to obtain a simple generalization of 
the canonical distribution in T-space by this approach. However, if we implement a mean- 
field approximation, a notion of generalized thermodynamics emerges in /z-space, as we 
have seen. In this context, generalized free energies are Lyapunov functionals associated 
with generalized Fokker-Planck equations. The fact that we cannot obtain a simple form 
of equilibrium distribution in T-space may suggest that the generalized thermodynamical 
formalism developed in /i-space is just effective. 



d_ 
<9r, ; 



D*pi 



. Pi 



C. Generalized Landau equation 

In [49], we have proposed to consider a generalized class of Landau equations of the form 

f = -( 2 -) D -^r / d^kk^umik ■ (v - voi//! (c»( f )§t - c "(fi)^)m 

where C is a convex function. For C = f In /, we recover the ordinary Landau equation 
(39) as a particular case. As shown in [49], Eq. (195) can be derived from the generalized 
Boltzmann equation introduced by Kaniadakis [46] in a weak deflexion limit. This equation 
conserves mass and energy and increases the generalized entropy S = — f C(f)d D rd D v. This 
corresponds to a microcanonical structure while the Fokker-Planck equation has a canonical 
structure as it decreases the free energy at fixed mass and temperature. 

In [49], we have also considered a test particle approach and showed the connection be- 
tween the generalized Landau equation and the generalized Fokker-Planck equation. Explicit 
expressions of the diffusion coefficient have been obtained in D = 3 for different entropies 
C(f). We consider here the case D — 1 and show that the generalized Fokker-Planck 
equation takes a relatively simple form. First, the generalized Landau equation in D = 1 
reads 

% = K± J dnS{v - VlWl (cV)%-CV^), (196) 
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where K is given by Eq. (47). Once again, we see that the collision term vanishes in D = 1. 
In the test particle approach, we must consider that f\ is fixed to the distribution of the 
bath, i.e. f\ = fo(v). This yields 



where D(v) = Kf (v). This equation governs the evolution of the density probability /(v, t) 
of a single test particle in a bath of field particles with fixed distribution fo{v). We note 
that, for t — > +00, the test particle relaxes to the distribution of the bath, i.e. / = f . 

V. CONCLUSION 

In this paper (Paper II), we have discussed the kinetic theory of Hamiltonian and Brown- 
ian systems with long-range interactions. The statistical equilibrium states of these systems 
have been considered in Paper I. For Hamiltonian systems, in the N — * +00 limit, we get the 
Vlasov equation. This equation admits an infinity of stationary solutions. One of them will 
be attained (on a coarse-grained scale) as a result of a violent relaxation on a short time scale 
of the order of the dynamical time ~ to- This metaequilibrium state, or quasi-stationary 
state (QSS), which depends on the initial conditions (due to the Casimir invariants) and 
on the efficiency of the mixing process (ergodicity), is usually difficult to predict [28]. The 
statistical equilibrium state, due to the development of correlations for finite iV systems, is 
reached on a much longer timescale t re i ax ~ N s tn increasing with the number of particles. 
For two and three-dimensional homogeneous systems, we can develop the kinetic theory 
at order 1/N and get the Lenard-Balescu equation. This equation converges towards the 
Maxwellian distribution establishing 5 — 1. For inhomogeneous gravitational systems, we get 
the Vlasov-Landau-Poisson equation when the collisions are treated as local and collective 
effects are neglected. This equation converges towards the mean-field Maxwell-Boltzmann 
distribution. Its dynamical stability (with respect to the Landau equation) coincides with 
the thermodynamical stability criterion in the microcanonical ensemble (entropy maximum). 
Because of logarithmic divergences of the diffusion coefficient, the relaxation time scales as 
Ueiax ~ ( N/ In N)t D ■ Finally, for one-dimensional systems (and 2D point vortices with 
monotonic angular velocity profile), the collision term vanishes at the order 1/N so that 
the evolution is due to higher order correlations. This implies 5 > 1 as observed by [12] 
for the HMF model. For Brownian systems, in the N — > +00 limit, we get the non-local 
Kramers and Smoluchowski equations. The mean-field Maxwell-Boltzmann distribution is 
the only stationary solution of these equations. Its dynamical stability (with respect to 
the Fokker-Planck equation) coincides with the thermodynamical stability criterion in the 
canonical ensemble (minimum of free energy). 

We have also formally introduced generalized kinetic equations and showed that they 
were associated with a generalized thermodynamical framework in \i space. We have intro- 
duced these equations (185), (187) and (192) from a specific class of stochastic processes 
(182)-(183) and (190) but we believe that these equations can have interest in much more 
general situations. They can be viewed as effective kinetic equations attempting to describe 
complex media. They may be useful when we are not in the strict conditions of applicability 
of standard kinetic theories. For example, there are situations in which the two-body dis- 
tribution function cannot be factorized as a product of two one-body distribution functions 




(197) 
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plus a small correction. In that case, the system is not described by the "ordinary" kinetic 
equations presented in Sees. II and III. On the other hand, the dynamics of complex sys- 
tems can be altered by microscopic constraints (hidden constraints) that modify the form 
of transition probabilities. Generalized kinetic equation can then be of interest to describe 
(at least heuristically) these non-ideal situations. One property of these generalized kinetic 
equations is to exhibit anomalous diffusion (since the diffusion coefficient depends on the 
density) and indeed, anomalous diffusion is observed in complex media. In that case, this 
is associated with a complex geometrical structure of phase space (e.g., multifractal in the 
case of the Tsallis entropy). Note that anomalous diffusion can also be due to the rapid 
decay of the diffusion coefficient with the velocity as in the Fokker-Planck equation (62). 
This is another, completely independent, reason for anomalous diffusion [11]. In that case, 
the Fokker-Planck equation is derived from the pure Hamiltonian dynamics and there is no 
relation with generalized thermodynamics. These two approaches apply to different regimes 
or different systems. 



APPENDIX A: THE LANDAU EQUATION AND THE TEMPORAL CORRE- 
LATION FUNCTION OF THE FORCE WITHOUT COLLECTIVE EFFECTS 

In this Appendix, we simplify the collision term appearing in the kinetic equation (38). 
For a homogeneous system, we need to compute the memory function 

dt J d D v 1 F ll (l -> 0, 0)F"(1 -> 0, t). (Al) 

The force created by particle 1 on particle at time t can be written 

F(l -> 0,f) = -im J ku(k)e ik(r(t) - ri(t Vk. (A2) 

Making a linear trajectory approximation r(t) = r + vt and ri(t) = ri + Vit, we get 

F(l -> 0,f) = -im J kw(k)e ik ( x+u Vk, (A3) 

where x = r — ri and u = v — Vi. Now, the memory function is easily calculated and yields 

AP" = n(2n) D m 2 J k"k v S(k ■ u)u(k) 2 d D k. (A4) 

Substituting this result in Eq. (38), we get the Landau equation (39). 



APPENDIX B: THE LENARD-BALESCU EQUATION AND THE TEMPORAL 
CORRELATION FUNCTION OF THE FORCE 

In this Appendix, we give a short derivation of the Lenard-Balescu equation. Its deriva- 
tion is classical but we shall need some intermediate steps in order to justify the expression 
(98) of the temporal correlation function. We follow the approach of Padmanabhan [3] but 
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we consider an arbitrary potential of interaction and we take into account collective effects. 
We start from the Klimontovich equation 



f + = »• (BD 

where f d (r, v, t) = J2i — i"j)5(v ~~ v i) * s the exact discrete distribution of particles and 
$d(r, t) is the potential that they generate. Writing f d — f + Sf and & d = $ + 5$ and 
taking the average of Eq. (Bl), we get for a homogeneous distribution 

| = £ ■ WV»>. (B2, 

Subtracting Eqs. (Bl) and (B2) and neglecting the quadratic terms (quasilinear approxi- 
mation), we obtain 

It can be shown that the terms neglected are of order 1/N 2 . We can solve Eq. (B3) by 
Laplace-Fourier transform. This yields 

df 

[u-k- v)Sf = -k • ^-5$. (B4) 
Therefore, the perturbed distribution function is given by 

5f(k, u, v) = &-6Q(Vi, u) + g(k, u, v), (B5) 

CO — K ■ V 

where g is the general solution of (co-k-w)g = 0. For g — 0, after integration over the velocity 
v, we obtain the dispersion relation e(k, u;) = which arises in the linear stability analysis 
of the Vlasov equation (see Sec. II C). The condition of marginal stability corresponds to 
e(k, 0) = 0. In the present context, g is related to the discrete nature of the system, i.e. the 
fact that the exact distribution function f d is a sum of ^-functions. Therefore, its expression 
is given by 

g(k, co, v) = — |p J2 m5(v - v,)e^(k • v - co). (B6) 

On the other hand, using the fact that the relation between the potential $ and the distri- 
bution function / is a convolution, we have in Fourier space 



5<Z>(k,oo) = (27r) D w(k) J 5f(k,oo,V)d D V. (B7) 
Combining Eqs. (B5) and (B7) we obtain 

<J*(k, co) = (2n) D ^- [ g(k, co, v)d D v, (B8) 
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where we have introduced the dielectric function (10). It is now possible to compute the 
collision current 

(5fV5$) = J (5/>,^,v)5<l(k^cJ0)ik , e ^(k+k '> r e- ^(wW)^ rf D krf D k , ^^ , . (B9) 
Using Eqs. (B5) and (B8), we get 



ik'{5f(k,u;,v)5$(k',u>')) = -(2tt) 



2D mu(k') 



e(k, uj)e{k!, u/) uj — k • v 



k' k 



dv 



e{k>,u>) 



From Eq. (B6), the correlation function is given by 

m 



x / (g(k,u;,v')g(k',u;',v"))d D Vd D V' + (27r) D — f^£(k')k' / (g(k,LU,v)g(k',Lu',v"))d D V 



(BIO) 



(g(k,u;,V)g(k',u; f ,v")) = ——f(V)5(v f - v")8(k ■ v' - cv)5(k + k')% + u'). (Bll) 



(27T) 

Substituting these expressions in Eq. (B9) and using 

P 



e(k,w) = l + (27r) D M(k) y 



a; — k ■ v 



m8{uj — k • v) 



av 



(B12) 



where P denotes the principal part, we finally obtain the Lenard-Balescu equation (49). 
Moreover, from the above expressions, it is easy to obtain the following expression for the 
correlations of the potential 



(5$(k, tv)5$(k', u')) = m(2n) D S(k + k')5(u> + u 



|e(k,c)| 2 

Taking the inverse Fourier transform for the time variable, we obtain 

<M>(k, 0)5*0*!, t)) = m(2n) D S(k + k>(k) 2 J |e(k;k \ v/)|2 
Noting that the force acting on the particle at time t is 

ik-vt jD\ 



f(V)5(k • v' - u)d D V(B13) 



ik - v,t f(V)d D V. (B14) 



F(t) = -J ik5$(k,t)e tk - vt d u k, 

we finally obtain the temporal correlation function of the force in the form 

u(k) 2 



(F^(0)F u (t)} = m(2n) u / k^k" 



|e(k, k • v 



-ik.(v-v') t/(v ^D v /^ k 



(B15) 



(B16) 



This expression generalizes Eq. (91) in the case where collective effects are properly ac- 
counted for. 
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